Maximizing genuine multipartite entanglement of N mixed qubits 
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Beyond the simplest case of bipartite qubits, the composite Hilbert space of multipartite systems 
is largely unexplored. In order to explore such systems, it is important to derive analytic expressions 
for parameters which characterize the system's state space. Two such parameters are the degree of 
genuine multipartite entanglement and the degree of mixedness of the system's state. We explore 
these two parameters for an iV-qubit system whose density matrix has an X form. We derive the class 
of states that has the maximum amount of genuine multipartite entanglement for a given amount 
of mixedness. We compare our results with the existing results for the N = 2 case. The critical 
amount of mixedness above which no iV-qubit X-state possesses genuine multipartite entanglement 
is derived. It is found that as N increases, states with higher mixedness can still be entangled. 

PACS numbers: 03.65.Ud, 03.67.-a, 03.67.Lx 



I. INTRODUCTION 

Quantification of entanglement of a bipartite state with 
arbitrary dimensions is a challenging problem and solu- 
tions are known only for some special cases with systems 
consisting of two qubits [l| , a qubit and a qutrit 0, Q , 
two qudits in a highly symmetric state [4L |5|, a double 
gaussian state jg], a pure bipartite state [3, Q, etc. Be- 
yond the bipartite case, quantifying the entanglement of 
a system consisting of three or more parties becomes even 
more challenging. For such systems, it becomes impor- 
tant to distinguish between genuinely multipartite entan- 
gled states, for which all the parties are entangled with 
every other party, from the bi-separable states for which 
some parties might be entangled with each other while 
being separable from the rest [9l-fllj|. 

In [l(| it was shown that a generalization of Wootters 
concurrence, which was derived for the two-qubit case 
[l[, is a good measure to quantify genuine multipartite 
entanglement. Although various witness operations to 
verify genuine multipartite entanglement existed fl2l - [l7j , 
it was not until recently that an analytic expression for 
Cgm was derived 18| . It was shown in 1181 that if a state 
of an A-qubit system has an X form one can get 
Cgm analytically. 

The A-states are density matrices whose only non-zero 
elements are the diagonal or the anti-diagonal elements 
when written in an orthonormal product basis of the 
A-qubit system. Although the properties of X density 
matrices for two-qubit systems are extensively studied 
in the literature [19l - l2l| . the properties of A-qubit A 
density matrices are still largely unexplored. In this re- 
port, we investigate the properties of these A-states of 
A-qubit systems. In general, an X-state of A-qubits has 
(2 +1 — 1) independent real parameters. Instead of ad- 
dressing this exponentially growing parameter space, we 
restrict our attention only to two fundamental quanti- 
ties, genuine multipartite concurrence and linear entropy. 
Linear entropy is a measure of mixedness and for a given 



state, say p, linear entropy is defined as 
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where d is the dimensionality of p. The linear entropy 
ranges from to 1 with indicating p to be a pure state 
and 1 indicating p to be completely mixed. In [22l |23| , 
Munro et al. explored the concurrence-entropy relations 
for all physically allowed two-qubit states. They found 
the class of states that has the maximum amount of en- 
tanglement for a given degree of entropy. It was shown 
that as the entropy increases, the maximum achievable 
amount of entanglement by any possible bipartite qubit 
state decreases. It was also shown in [22l - |24j that beyond 
a critical value of S cr — 8/9, no two-qubit state could 
possibly be entangled. The class of states that have the 
maximum amount of entanglement for a given degree of 
entropy is referred to as maximally entangled mixed states 
(MEMS) [11]. For the two-qubit case, the MEMS class 
happens to be of the A form. This fact suggests that 
even for an A-qubit system, they may be of the A form. 

Here, restricting our analysis to the class of A-states, 
we will generalize the two-qubit results to a general N- 
qubit setup. In particular, we find the class of A density 
matrices that, for a given amount of entropy, has the 
maximum amount of genuine A-qubit entanglement (see 
Eq. (flfj)) ). We find that the critical amount of entropy, 
as a function of A, above which no A density matrix can 
possibly exhibit genuine multipartite entanglement is 



12JV-1 



Scr(N) 



(2 N - 1)(2 N - 1 + 1)' 



(2) 



We see that this critical entropy increases as a function 
of A. This indicates that in the limit of a large num- 
ber of qubits, it is possible to find highly mixed states 
that can still possess entanglement. The dependence of 
various physical parameters on the entropy of a bipartite 
system state and the scaling of these dependences on the 
dimensionality of the two parties involved can be found 
in MM. 
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In Section UH we briefly explain the concept of gen- 
uine multipartite entanglement. In this section, we also 
describe the A density matrices of A qubits and the gen- 
uine multipartite concurrence that quantifies the entan- 
glement for these states. In Section IIII1 we derive the 
maximally entangled mixed A states (X-MEMS) for N 
qubits. We conclude in Section Hvl 



II. GENUINE MULTIPARTITE 
CONCURRENCE OF iV-QUBIT X DENSITY 
MATRICES 



A system consisting of TV qubits is said to possess gen- 
uine multipartite entanglement if each qubit is entangled 
with every other qubit and not only to some of them 
To make this statement more explicit, we first 
define what is meant by a bi-separable state: A pure 
A-qubit state, say \ip), is said to be bi-separable if it 
can be written as a product of two multi-qubit states: 
\i>ba) — \ipi)®\*p2) (the subscript stands for bi-separable) . 
A mixed state is bi-separable if it can be written as a mix- 
ture of bi-separable pure states: pb s = SjPi|V , L)(V'Ll- 
If an A-qubit state is not bi-separable, it is genuinely 
iV-partite entangled [91-fTlj. 

One of the proposed measures to quantify genuine 
multi-partite entanglement is genuine multi-partite con- 
currence [nj. For a pure state, genuine multi-partite 
concurrence is defined as: 



C GM (m := mm ^2(1 - Tr(p* )), (3) 
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where Y represents the set of all possible bi-partitions, 
{A 7 |i? 7 }, of the N parties and pi is the reduced density 
matrix: p& — Tr# T (\ip)(ip\). For a mixed state, p, the 
concept of Cqm can be formally generalized to 



eral structure of the following form: 



Cgm(p) = inf 

{p*,\r)} 



(4) 



where the minimization has to be carried out over all 
possible pure state decompositions of the density matrix: 
p = ^liPi^) {il> % \- This minimization procedure renders 
analytic or even numerical parameterization of Cqm in- 
feasible in general. 

For some cases, lower bounds on Cqm are available 
for analysis [13 - [l7| . Only recently, by generalizing the 
result for the two-qubit case [l| , an exact analytic expres- 
sion for Cqm was derived for an A-qubit system whose 
density matrix has an X form 18j. To briefly elaborate 
on the result in [l8j], let us look at an A-qubit X den- 
sity matrix. Arranging the rows and columns in an or- 
thonormal product basis, |1, 1, • • • , 1, 1), |1, 1, • • • , 1, 0), 
• • • , |0, 0, • • • ,0,0), the X density matrices have the gen- 
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where n — 2 A '~ 1 . For the X matrix to be a valid density 
matrix, one must have X)i( a i + M = 1 an d \zi\ < VaJh- 
These conditions are necessary and sufficient for the X 
matrix to have unit trace and positive eigenvalues. In 
18], it was shown that the Cqm for the above A-qubit 
state is given by 



Cgm{X) 
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Note that for the two-qubit case, this result reduces to 
the previously known result derived by Wootters [l[ . 

If any one of the qubit's state is traced out from the 
A-matrix, the rest of the qubits become separable. This 
fact can be understood by noticing that if one traces 
any of the qubits from Eq. (0, one gets a diagonal and 
hence separable reduced density matrix. An example of 
A-states is the class of well known GHZ states. The GHZ 
states have theoretical and practical importance. In [28[ , 
it was shown that there are no bounds to the amount 
by which the GHZ states can violate the limits imposed 
by a Bell's inequality. In the context of trapped ions, it 
has been possible to experimentally create GHZ states of 
up to fourteen qubits [29(. With the advances in experi- 
mental circuit QED, it is predicted that large number of 
qubits can be initialized in the GHZ states |3Cf . For two- 
qubit systems, the A-states are central to understanding 
the relations between fully entangled fractions and en- 
tanglement [3l| . Taking into account the experimentally 
important scenario of un-controlable environmental in- 
fluences on the system, it was shown in [3, [l!| that the 
A-state form is robust in the sense that states initially 
in the A form remain in the X form if local damping 
channels act on the various non-interacting qubits. The 
possibility of an arbitrary initial pure state de-cohering 
into an A-state was studied in [20]. Thus, we see that 
the A-states are not just a mere mathematical construc- 
tion but are an important class of all physically allowed 
A-qubit states. 



III. MAXIMALLY ENTANGLED MIXED 
A-STATES FOR iV-QUBITS 

In this section, we explore the entropy verses gen- 
uine multipartite entanglement properties of A-qubit A- 
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states. For concreteness, we begin by considering a three- 
qubit system (The two-qubit system was studied in f22j). 
We randomly generate a million X-density matrices of 
three qubits and determine their linear entropy and gen- 
uine multi-partite concurrence using Eqs. ([T]) and (O 
respectively. The entanglement-entropy plane for the 
three-qubit X-states is shown in Fig. [T] where we have 
plotted Cqm and S for each of the randomly generated 
X matrices. We notice two qualitative trends from Fig. 



the matrix elements of X are as follows: 
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FIG. 1. (Color online) Genuine multipartite entanglement 
and entropy plane for three-qubit X-states. Each of the mil- 
lion points corresponds to the Cgm(X) and S(X) calculated 
for a randomly generated X-density matrix. We are inter- 
ested in the class of states that lies at the boundary. We also 
want to know what the critical entropy, S cr , is as a function 
of N. Critical entropy is the entropy above which no possi- 
ble system state in the X form possesses genuine multipartite 
entanglement. 



[TJ First is that as the entropy increases, the maximum 
amount of achievable genuine multi-partite entanglement 
decreases. The second important thing to notice is that 
beyond a critical amount of entropy (close to 0.9 in this 
case), no state possesses genuine multipartite entangle- 
ment. These qualitative features can be quantitatively 
captured by knowing the class of states that lies at the 
boundary of the entanglement-entropy plane. These are 
the states that have the maximum entropy for a given 
amount of entanglement or said in a different way, has the 
maximum amount of entanglement for a given amount 
of entropy. We will refer to this class of states as maxi- 
mally entangled mixed X -states (X-MEMS). The purpose 
of this report is to find X-MEMS for an TV-qubit system. 



Without loss of generality, we assume that for a gen- 
uinely multi-partite entangled state the concurrence is 
given by, C GM {X) = 2(\z 1 \ - YTj=t\ \/ a ih)- We now 
prove that for every density matrix of the form given in 
Eq. ([5]), we can find another density matrix, X, that has 
the same Cqm as the density matrix X but has a higher 
entropy. The density matrix elements of X in terms of 
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The 6i's are added to the a^'s to ensure that the trace 
of X is one. It is clear that the genuine multi-partite 
entanglement of X and X is the same: 



Cgm{X) — Cgm(X). 



(8) 



In order to prove that the entropy of X is more than the 
entropy of X, we note that for \zi\ > Y^i we have 



S(X)-S(X) 
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The facts that C GM (X) = C GM {X) and S{X) > S(X) 
imply that the maximally entangled mixed X states have 
to be of the X form for which all the lower diagonal 
elements except b\ are zero. 

To reduce mathematical complications while deriving 
the form of X-MEMS, we redefine the following matrix 
elements of X which were previously defined in Eq. ([7]): 



ai = a + |7|, 

h = /3 + h\, 
zi = 7. 



(10) 



To ensure the positivity and unit trace properties of the 
density matrix X, the following conditions must be sat- 
isfied: 

a + W >0, 

P + >o, 

ap + ia + PM >0, 



a + /3 + 2|7| = 1. 



(11) 



i=2 



With this redefinition, the density matrix, X, takes the 
following form: 

/a + \j\ 

a-2 



X = 







7 \ 
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V 7* 



(12) 



4 



For the above form of the density matrix, the genuine 
multipartite concurrence is simply given by Cgm(X) = 
2\j\. Keeping [7] fixed, we now need to figure out the 
expressions for a.;, a and (3 that would maximize the en- 
tropy, S(X). To this end, we note that because of the 
presence of the off-diagonal element 7, the diagonal ma- 
trix elements a and j3 have to be treated separately than 
all other diagonal entries, ai for 2 < i < n. We also note 
that keeping the off-diagonal element constant, entropy 
or mixedness is maximized if all the possible diagonal 
elements are equally populated, i.e. when 
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y, 



(13) 



where we have defined two new variables, x and y. Using 
the constraint that the sum of the diagonal terms should 
sum to one, the variable y can be eliminated from the 
expression of the linear entropy of the density matrix X 
to get: 



S(X) 



where 



2 N _ 



T (A(x+| 7 |) 2 + S(x+|7|) + c), (14) 



A= -2 



D 



n + 1 
n-l) ' 



n- 1 



C= 1- 



n — 1 



-2| 7 | : 



(15) 



We now need to maximize the linear entropy, S(X), for 
a given amount of Cgm(X) — 2\j\. While doing so, we 
have to remember the constraints in Eq. (1111) . According 
to Eq. (|TT|) . one must have x,y > and 2(x+ | 7 |) + (n — 
l)y = 1, which implies that (x + | 7 |) < 1/2. Keeping 
these constraints in mind, we plot S(X) in Fig. [5] as a 
function of x + | 7 | for four different values of | 7 |. We 
note that when I7I < l/(n + 1) the maximum of S(X) 
occurs at [x + I7I) = l/(n + 1), which according to the 
trace condition means that y = l/(n+ 1). On the other 
hand, when I7I > l/(n+l), the maximum of S(X) occurs 
at (x + I7I) = I7I, which in turn implies that y = (1 — 
2|7|)/(n — 1). With these considerations, we finally get 
the following form of the A-density matrices that has 
the maximum amount of entropy for a given amount of 
genuine multipartite entanglement: 



//(7) 
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FIG. 2. (Color online) Maximizing the entropy, S(X). The 
solid lines correspond to S(X) for different values of j — y | over 
the allowed range of (x + According to Eqs. and 
(fT3)) . the allowed range of (x + is | 1 < (x + | 7 |) < 1/2. 
The solid dots correspond to the respective maximas of the 
various curves. We see that the maximum of S(X) occurs at 
(x + | 7 |) = l/(n + 1) if | 7 | < l/(n + 1) and at (x + | 7 |) = | 7 | 
if | 7 ] > l/(n+l). 



where 



l/(n+l) < | 7 | < l/(n+ 1) 
l7l V(n + 1) < l7l < 1/2, 



(17) 



and 
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l/(n + l) 
(l-2| 7 |)/(n-l) 



< I7I < l/(n + l) 
l/(n+l) < | 7 | < 1/2. 

(18) 



As noted earlier, these are also the states that have the 
maximum amount of entanglement for a given amount 
of entropy. The degree of genuine multipartite entangle- 
ment corresponding to these states is Cgm(X) = 2\j\ 
and the corresponding linear entropy is: 



S(X) - w — i (^/ 2 ( 7 ) + B/( 7 ) + C) 



(19) 



The form of the matrix X derived in Eq. (|16p in terms of 
the concurrence, Cgm — 2| 7 |, is the main result of this 
report. This class of states is the maximally entangled 
mixed A"-states (X-MEMS). To illustrate our analytically 
derived results, we plot in Fig. [3] the entanglement ver- 
sus entropy of a million randomly generated three-qubit 
and five-qubit A"-density matrices. The analytical bound 
found in Eq. (fTfjj) is plotted as a solid line. In clear agree- 
ment with the analytic bound imposed in Eq. (|19l) . we 
see that all the randomly generated points lie inside the 
solid line. 

It should be noted that the general results derived here 
for the -ZV-qubit case reduces correctly to the results de- 
rived by Munro et al. for N = 2 [221] . We reiterate the 
fact that the analysis in this report was restricted to a 
special form of iV-qubit density matrices that have the 
X form as given in Eq. ([S]). On the other hand, the re- 
sults of [22| were valid for all physically allowed two-qubit 
states and not just for A-matrices. 
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FIG. 3. (Color online) Entanglement and entropy of randomly 
generated three-qubit and five-qubit X-density matrices. The 
solid lines correspond to Eq. (|19p which analytically gives the 
maximum entropy for a given amount concurrence. We see 
that in agreement with the analytical results derived in this 
report, all the random points lie inside the solid line. Note 
that as the number of qubits increases from N = 3 to N — 5, 
the critical amount of entanglement beyond which no state is 
entangled also increases. 



Putting I7I = in the Eq. (|19|). we can get the critical 
amount of entropy beyond which no A-states can exhibit 
entanglement. This critical entropy as a function of A 
turns out to be: 
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which for N = 2 correctly gives S cr (2) = 8/9 [22j. In Fig. 
21 we plot S cr (N) as a function of the number of partic- 
ipating qubits, N. We see that S cr (N) is an increasing 
function of A and approaches unity asymptotically. This 
implies that as the number of qubits increases, we can 
have states that are more mixed and still exhibit genuine 
multipartite entanglement. 



IV. CONCLUSION 

In this report, we found the class of A-qubit A-states 
that has the maximum amount of genuine multipartite 
entanglement for a given amount of mixedness. This 
class of states derived here, referred to as maximally en- 
tangled mixed X-states (X-MEMS), is an extension to 
the ones derived for the two-qubit case in (22|. We find 
that the critical amount of entropy, S cr (N), above which 
no X-state possesses genuine multipartite entanglement 



increases as a function of the number of participating 
1i — 
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FIG. 4. (Color online) Critical entropy, S cr (N), as a function 
of the number of participating qubits, N. We see that the 
critical entropy increases as the number of qubits increases 
and asymptotically approaches unity. 



qubits, A. As the derived class of states has the maxi- 
mum amount of entanglement for a given amount of en- 
tropy, they might have practical utility in cases that re- 
quire entanglement as a resource but are influenced by 
unavoidable sources of noise. 

The entire analysis in this report was restricted only to 
density matrices that have an X form. This restriction 
was primarily because of the fact that a computable mea- 
sure which captures genuine multipartite entanglement 
of an A-qubit system exists only for X-density matrices. 
The results in [22| for the two qubit case were more gen- 
eral. In [22|], maximally entangled mixed states (MEMS) 
were found out of all possible two-qubit density matrices 
and the analysis was not just restricted to density matri- 
ces of the X form. The MEMS found in [22| happens to 
be of the X form. This fact suggests that the MEMS of 
an A-qubit system might also be of the X-form. If this 
happens to be true, the MEMS of an A-qubit system 
will be the same as the X-MEMS found in this report. 
Since an analytic form of Cqm for an arbitrary A-qubit 
density matrix is yet to emerge, a connection between 
MEMS and X-MEMS of an A-qubit system remains an 
open question. 

For the two-qubit case, it was shown that the maxi- 
mally entangled mixed states depend upon the measures 
used to quantify mixedness and entanglement. In this 
report, the measure of mixedness was chosen to be linear 
entropy (see Eq. (JTJ), and the measure chosen to quan- 
tify entanglement was genuine multipartite concurrence 
(see Eq. ©). Since Cqm is the only computable measure 
available to quantify A-qubit entanglement, the depen- 
dence of X-MEMS on different measures of entanglement 
is still an open question. 
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